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Abstract. - The measured and reported rate [E| of the 3 He A — > B transition 
depends on the diameter of the experimental cell. Theoretically, this can only be 
explained by a direct drag on the interface exerted by the lateral walls, leading to a 
pressure discontinuity across the interface. This has not been considered before, ei- 
ther here or in similar circumstances. The general hydrodynamic theory that includes 
the influence of finite geometries is presented. 
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Introduction. - General hydro dynamic considerations |2|-|J on the 3 He A — > B transi- 
tion have in the past yielded a fairly complete picture for the infinite geometry. Depending 
on the interface velocity u, there are two qualitatively different types of behavior. If u is 
small compared to the second-sound velocity, « < C2, three characteristic features of the 
phase coherent transition are predicted: (i) The chemical potential is continuous across the 
interface, forcing the melting A phase to be warmer than the growing B phase, (ii) The 
relevant mode of heat transfer, second sound, is efficient enough to obliterate any onset of 
hypercooling. (iii) The effective growth resistance preventing the transition rate to diverge 
consists of three elements in series: two account for the sq-mode, a hydrodynamically gentle 
temperature variation on each side of the interface, the third accounts for a sharp tem- 
perature discontinuity at the microscopically narrow interface. The former two dissipative 
elements stem from collisions among quasiparticles, the latter one (probably much smaller 
in magnitude) results from Andreev scattering of these particles at the interface. At the fast 
end of the transition, ii > C2, circumstances are rather different: (i) Second sound is now 
a comparatively slow process, leaving the chemical potential discontinuous and the growing 
B phase warmer, since an ever-increasing portion of the latent heat is no longer removed 
from the interface region, (ii) The usual scenario of hypercooling (in which no heat at all 
is transfered) becomes asymptotically valid, (iii) The total growth resistance contains both 
additive and multiplicative contributions. 

Since second sound is so important in determining the properties of the A-B transition, 
it seems plausible that the above features would not easily survive in a confining geometry: 
walls hold back the normal velocity v n and thus confine the range of heat transfer via second 
sound. In what follows we briefly state what this range is, why its finiteness is important, 
and that it alone is insufficient to account for the measured |]J geometry dependence of the 
phase transition. Generally, for stationary flows, lateral walls are well-accounted for by a 
relaxation term pv n /r in the Navier-Stokes equation (obtained by averaging over the cross 
section of the experimental cell, eg. r = pR 2 /8rj for a cylinder of radius R). As a result of 
this term, the second-sound step functions sent out by the moving interface p|-§J acquire 
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an extra, exponentially decaying factor. In the rest frame of the interface, with B phase at 
x < and A phase at x > 0, it is exp(x/A2). The resulting stationary temperature profile is 
shown in Fig. 1 and 2. (On the scale of these figures, the temperature variations of the sq- 
and diffusive modes [|J are too small to be separated from the discontinuity at the interface, 
x = 0.) The decay length 

A 2 = -p n r(cl - u 2 )/(p s u) (1) 

is also the range of heat transfer via second sound. Outside this range, |x| 3> [ A2 [ , there is 
no means for heat transfer, and energy balances locally. Now, it is important to realize that 
the usual thermodynamic consideration of hypercooling is based on local balance of energy. 
Hence, the temperature T B = T A + 5T of the growing B phase is taken to be given by the 
A phase temperature T A and the heating 5T from the released latent heat. Second sound 
transfers heat away from the interface and invalidates this consideration. However, if the 
heat transfer has a finite range, hypercooling remains a valid concept outside this range. 
And indeed, Fig. 1 shows that the B phase is warmer than the A phase for x 3> [ A2 1 - In 
addition, the solution depicted in Fig. 1 is obtained under the condition 

'\Afi\ - (Aa) 2 /d T a 

which is, up to first order in Tab—T, equivalent to the usual and approximative hypercooling 
condition (Tab ~ T) — Iab/c > 0. (T : initial undercooling temperature; Tab '■ coexistence 
temperature; \i : chemical potential, a : entropy, Iab '■ latent heat, c : specific heat, all per 
unit mass; A : difference across the interface, eg. A/i = ji B — fi A .) 

Above the hypercooled region, e = 1 — T/Tab ~ 0.5%, planar instability and non- 
stationarity (the usual conundrum of snow flakes) may be expected to be relevant, but are 
not: Approaching the critical value e c ~ 0.5 % from below, u vanishes with the left side of 
Eq. (0), and A2 diverges accordingly. Once A2 exceeds the longitudinal dimensions of the 
experimental cell, even a heat transfer within the scale of A2 is sufficient to disturb the energy 
balance and invalidate the concept of hypercooling. Then a planar transition becomes viable 
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again. (The appropriate solution is non-stationary, similar to that of Ref . ; see discussion 
at the end of the introduction.) 

Concentrating now on the hypercooled region, where the stationary solutions of Fig. 1 
and 2 are valid, we nevertheless face the following difficulty: Despite the relaxation term 
pv n /r and the resulting qualitatively changed behavior of second sound, the interface veloc- 
ity u remains independent of the cell radius R, in stark contrast with the well-documented 
experimental observations Q. This is indeed perplexing, since the hydrodynamic descrip- 
tion accounts for the structure of the dynamics and is independent of specific mechanisms, 
especially if any has been overlooked. Re-examining the equations of Ref. we found 

no reason to doubt the strict conservation of energy and mass, leaving the momentum as 
the only candidate to resolve the puzzle: Rough lateral walls break the translational sym- 
metry and it is conceivable that the Navier-Stokes equation acquires a relaxation term that 
is singular at the interface, resulting in a discontinuity of the momentum current, 

Ap = Tpu/R . (3) 

This equation can be viewed as the equation of motion for the interface, in which the 
inertial term ~ u and the restoring force ~ u vanish; the pressure discontinuity Ap is the 
external force (per unit area), while the damping ~ ii comes from the drag exerted by the 
lateral walls and would have its origin, for instance, in the interaction between the surface 
tension (including textural contributions) and the wall's roughness. The total drag should 
be proportional to the contact length, 2irR in a cylinder, per unit area it is ~ 1/R. V is 
introduced as a phenomenological parameter, accounting for the degree of roughness. As 
we shall see, Eq. (|3]) indeed establishes agreement between theory and experiment. More 
specifically, we find a "wall contribution" ~ 1/R added to the total growth resistance of 
the infinite geometry, so that the change of ii is also proportional to 1/R. (The necessity of 
having Ap ^ has not arisen before, such as at the well-studied 4 He II-vapor or 4 He Il-solid 
interfaces |5|]. At least, it was not studied there. A re-examination may be worthwhile.) 
Finally, we address the relationship between the stationary solution of second sound 
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presented here and the propagating step functions of Ref. in the case of an infinite 

geometry. It is noteworthy that letting R^oo, the exponentially decaying solution does not 
reduce to the same steps. This is because of the following reasons: In deriving the stationary 
solution here, the lateral dimension is taken to be finite, the longitudinal one infinite; and 
although an arbitrarily large R results in an arbitrarily large decay length A2, it is finite 
and always smaller than the longitudinal dimension. The energy current Q in the present 
solution is constant, it always is in stationary, one-dimensional solutions, d x Q = —dte = 0. 
Since Q vanishes outside A2 (aside from a constant equilibrium contribution) it must also 
vanish within |J. This alone is enough to exclude the type of solution represented by 
the step functions propagating toward infinity, all the while transferring energy. The two 
types of solutions are in fact mutually exclusive: If the lateral geometry is confining enough 
to entail a decay length A2 much smaller than the longitudinal dimension, the stationary 
solution presented here is relevant, while the propagating solution is valid in the opposite, 
open limit, where the A2 is much larger than the longitudinal dimension. 

Theory. - In what follows, we shall derive the general stationary solution of superfluid 
hydrodynamics for the temperature T, the counterflow velocity w = p s (v n — v s )/p n and the 
A-B interface velocity it, taking into account the relaxation term pv n /r in the momentum 
conservation, and the altered connecting condition, Eq. (H). The solutions for T and w we 
search for are stationary in the rest frame of the interface, with the B phase at x < and 
the A phase at x > 0. Before the actual calculation, however, we need to point out the 
following consideration that significantly simplifies our task: As long as the radius of the 
experimental cell is not too small, R > 0.1cm, the physics close to the interface, especially 
the sq- and the diffusive lengths, are little changed from the infinite geometry, and most of 
the results from Ref. 0j remain valid and need not be recalculated. More specifically, we have 
three characteristic length scales: (i) The microscopic width of the interface, too small to be 
considered in a hydrodynamic approach except as the location of discontinuities, (ii) The 
large decay length of second sound, Eq. ([I]), A 2 ~ — A/jc 2 /-u for «i<c 2 and A 2 ~ \ru/c2 for 
ii ^> c 2 , where the characteristic length A# = rc 2 p n / p s is of order 0.1cm for radii R ~ 0.1cm. 
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(iii) The mesoscopic decay lengths of various exponential functions. Although individually 
describable in a hydrodynamic theory, they may also be collectively accounted for by a total 
growth resistance of an effective interface |4]] (that is wide enough to include these decays). 
The decay lengths are the sq-length X sq for ?i C C2 and the diffusive lengths Ay, X® for 
u ^> C2. In an infinite geometry, all three are combinations of two characteristic lengths, both 
of order KTW A t = k/(2c 2 pTd T a) and X w = [(4/3) tj - p(& + CO + C2 + p 2 Qps/ (2pp n c 2 ). 
(The heat conductance k and the viscosities 77, C1-4 are defined in the usual way [0], neglecting 
anisotropy.) In a finite geometry, the mesoscopic decay lengths are modified in first order of 
At/ Xr, X w / Xr as 

and 

X? = 2Ax~r f 1 - t^V^4) , X S = 2A,^ f 1 4- ^VM) • (5) 



The correction terms obviously may be safely neglected for radii R ^ 0.1cm. For this ex- 
perimentally relevant case (the only one we study M), we may therefore take the effective 
interface as calculated in Ref. [§]] for the infinite geometry, especially the total growth resis- 
tance, and only consider the change further away from the interface, on length scales of Xr. 
(Technically, this allows us to discard all dissipative terms in the hydrodynamic equations 
except the new relaxation term pv n /r, cf. the detailed discussion on the effective interface 
in Ref. ||.) 

Linearizing the averaged hydrodynamic equations with respect to w = p s (v n — v s )/p and 
T A ' B — Tj, and retaining terms of second order in u/c 2 an d 02/11 respectively, the stationary 
solution for u <C c 2 is given by 

T B (x) = Tj + 5T B , T A (x) = T t + ST A exJ--^-) , (6a) 

u x 

C2 X R/ 



w B (x) = , w A (x) = 5w A expf--^- ) , (6b) 



and for ix ^> c 2 by 



T B {x) =T i + 5T B + ST B exp(^) , T A (x) = T { , (7a) 

V u Any 

W B ( X ) = 5w B exp( -^-) , w A (x) = . (7b) 
\u\rJ 

As discussed above, second-sound and the related energy transfer now acquires a finite 
range. Its decay length A 2 = — Xr(cI — ■u 2 )/(c 2 , u) is negative for u < c 2 and positive for 
u > c 2 . (Therefore, ST B , Sw B = in the first and ST A , 5w A = in the second case.) In a 
confined geometry, the amplitudes of the stationary second-sound modes are related by 

5w A ' B = u (drcr/cr) 5T A ' B (8) 

(in contrast to the well-known relation 5w A,B = ±c 2 (dxcr/a) 5T A ' B for the infinite geome- 
try). The constant 5T B accounts for the temperature rise due to the latent heat and the 
interface dissipation, while 5T A = due to causality and the lack of any means for long 
range energy transfer. The temperature fields (outside the effective interface), are plotted 
in Fig. 1. and 2. 

The interface velocity u and the independent amplitudes of Eqs. ([]), (0) are to be de- 
termined from the general connecting conditions (CCs) at the effective interface: 



AQ = , A# = , Ap = Ypii/R , (9) 
A(ji + v n v.) = Q, Rs = (f) AT + gAfi. (10) 

All quantities are defined in the rest frame of the interface and have been averaged over 
the cross section of the experimental cell. The mass current is given as g = —pit @3, (/) 



denotes the average \{f B + f A ) of the entropy current across the interface. Except the 
non- vanishing discontinuity of the pressure p, all CCs are as derived in Ref . || , though the 
subscript e for "effective" has been eliminated from A and ( ). From the derivation there, 
it is also clear that R s is the total entropy production within the effective interface. The 
entropy produced by the drag exerted by the lateral walls at the rate of iiAp = Tpu 2 /R is 
therefore implicitly included in R s . Subtracting this part, one obtains the rate of entropy 
produced in an infinite geometry 
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Rf — R s — Tpu 2 /R , (11) 

which was parametrized as Rf* = (f) 2 /n e and g 2 /K e , respectively for ii <C c 2 and it ^> c 2 ; 
cf. [fj|. (Explicit expressions for the effective, R-independent Onsager coefficients K e , K e can 
be found in Ref. [f|.) Since the term Tpit 2 /R in Eq. ([TT|) can be rewritten as T(f) 2 / (Rp(a) 2 ) 
and Tg 2 /(Rp), we finally obtain from R s ~ (/) 2 and g 2 the CC: 

(/) = (^ 1 + r/( J Rp( ( x) 2 ))^AT 



and 



g = (K- 1 + T/(Rp)Y (Ap Q + (a) AT) 



(12) 
(13) 

respectively for u <C c 2 and u 3> c 2 . Here, p a denotes the chemical potential for given 
pressure and temperature, in a system with v n — v s — 0. 

Expanding the CCs to first order in the deviation from the initial temperature Tj, we 
obtain a set of equations for the interface velocity it and the amplitudes 5T B , 5T A (or Sw B ). 
In the limit of ii <C c 2 , it is solved by 
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(14a) 
(14b) 

(14c) 



while the dependent counterflow amplitude is given as 5w A = u[dT<y A /<y A ]iST A . Note 
the independence of the constant temperature 5T B , caused by the drag ~ T. Since u is 
taken to be positive here (melting the A phase on the right), the validity of Eqs. (|T4| ) is 
confined to initial temperatures satisfying [|A/i | — (Aa) 2 / 'dr(y B }i > 0, cf. the discussion in 
the introduction. For u 3> c 2 , expansion of the CCs around Tj and linearization with respect 
to c 2 /u, ST B and 5w B lead to 

Act Apo + TuR- 1 ' 



5T 1 



Sw B 



u 



d T a B Td T a B 
p s Ap + TiiR- 1 - a B 5T B 



Pn 

p 2T 



u 



R 



-1 



1 



Ap - -Ao5T B 



(15a) 
(15b) 
(15c) 



while the dependent temperature amplitude is given as 5T B = {l/u)[a B / dT(T B ]i Sw B . Again, 
5T B contains an .R-dependent contribution from the latent heat and the interface entropy- 
production R s ; cf. Eq. (|14]a). 

Summary. - The influence of finite geometries is considered, the interface velocity u 
and the stationary temperature field are derived. For both u <C c<i and u 3> C2, the total 
growth resistance of the transition is a sum of the open geometry Onsager-coefficient and 
a "wall contribution" ~ T/R. Hence u (for given Tj) is maximal for R — > 00, as measured 
Q. The phenomenological parameter T depends on the cell's inner surface, cf. Eq. ([|). 
Two sets of it, obtained in tubes of different radii R, should suffice to determine T. This 
is not possible with the published data 0, however, because the tubes used were made of 
different materials, copper and epoxy, with presumably rather different T. Perhaps further 
experimental will fill this gap. 

Defining a total growth coefficient as K tot {Ti/TAB , H) = (l/K e + 2T/Rp)~ 1 , we esti- 
mate from the data of Ref. [| that K tot (0M , 1.5kOe)= 1.9gscm" 4 , K tot (0M , 3.0kOe) = 
2.8gscm" 4 , K tot (0.85, 1.5kOe)= 1.3gscm- 4 , K tot (0.85 , 3.0kOe)= 1.4gscm- 4 ; obvi- 
ously, K tot increases with increasing magnetic field or decreasing supercooling rate. 

A paper containing the algebraic details of this letter and of Ref. || is being prepared. 
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FIG. 1. The temperature field for u <C c 2 , as in Eq. (6a). 
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FIG. 2. The temperature field for u 3> c 2 , as in Eq. (7a). 
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